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OBJETIVOS 


1. Entender y resolver problemas relacionados con problemas de 
conduccion de calor en regimen transitorio 


2. Aprender a resolver ecuaciones en derivadas parciales por 
medio de metodos numericos. Aplicadas a la resolucion de 
fenomenos de conduccion (difusion) en regimen permanente y 
transitorio. 
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MANUALES DE PROGRAMACION 


1 . http://www.programatium.eom/c.htm 

2. http://www.cplusplus.com/doc/tutorial/ 

3 . 

http://www.cprogramming.com/tutorial.html 

4 . Wikipedia... 
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CONDUCCION TRAN SITORI A 


■ En el curso se ha ido incrementando gradualmente la 

complejidad de las situaciones estudiadas. Pero no se han 
estudiado casos en los que la temperatura varie con el tiempo 
(transitorio) 

• Los problemas de transferencia de calor dependientes del 
tiempo son frecuentes. Surgen debido a la variacion de las 
condiciones de frontera 

■ Al inicio se considerara que la temperatura dentro del solido 

permanece constante T(t), luego la temperatura podra variar 
con la posicion y con el tiempo T(x,t) en 1 D o T(x,y,t) en 2D 
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Sist. concentrados/resistencia interna despreciable 


Los sistemas concentrados son aquellos en los que la 
temperatura interior no varia con la posicion: T(t)==> todo el 
sistema se encuentra a la misma temperatura 


* Bola de cobre vs trozo de carne 


Ocurre cuando la resistencia al interior del cuerpo es muy 
pequena comparada con la resistencia a la transferencia de 
calor entre el sistema y los alrededores. 


• La transferencia de calor se determina por medio de un 
balance de calor global 

( Transf. calor desde el cuerpo ) = ( Increm. energia en el tiempo ) 


-hA(T-TJ] = mc dT 


oo 


dt 
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Reemplazando m= p V y haciendo un cambio de variable 
6 = T-T 


oo 


-hAO = pV c 


d_0 

dt 


e 


-pVc r de. 


j^=j* 


pVcde , 

m~T = * 

_-p Vc 


hA q 0 0 


t = 


hA 


in(e/e,.) 


O tambien 


e 

e. 


t 


/iA 


= e p Vc 


El inverso de la cte en el exp 
se conoce como cte de tiempo 
termica Tt 
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Conociendo la temperatura se puede evaluar el intercambio de calor 
con los alrededores y por tanto tambien el cambio en la energia 
interna del sistema 

Q{t)=hA{T{t)-Tj [W] 

La cantidad total de calor transferida entre el sistema y el medio 
desde t=0 hasta t=t sera 

Q{t) = mc{T{t)-T I ) [kJ] 

El calor maximo transferido se alcanzara cuando la temperatura del 

cuerpo sea la de los alrededores 

Qma X = m c{T x -T I ) [kJ] 
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* Dada la sencillez del metodo resulta importante establecer 
cuando es adecuado aplicarlo 


• Existe un parametro que cuantifica cuan exacto es el metodo. 
Para obtenerlo se considera la conduccion en una placa que es 
enfriada por conveccion. 


t 



T —T 

J<A - — — s ——hA{T s2 —T 

J—i 


oo 


T s ,-Tj_(L/kA)_R 


cond 




00 


1/M R 


conv 


T s- 1 Ts ’ 2 ’=hk= Bi 
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R 

R 


cond 


conv 



Bi es el numero de Biot. Por conveniencia 
se define L=V/As (facil geom complejas) 


Puede ser interpretado como la relacion 
entre las resistencias. Si Bi« 1 la resistencia 
de conduccion es mucho menor que la de 
conveccion. El metodo puede ser aplicado. 


r t 



Tambien se puede entender como la relacion entre la diferencia de 
temperaturas en el solido, respecto a la de los alrededores 


Si Bi« 1 la variacion de la temperatura dentro del cuerpo es pequena 
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En general se considera adecuado el uso del metodo de resistencia 
despreciable o modelo de sistema concentrado si 

Bi<0.1 

Lc: distancia de la maxima 
diferencia de temperatura 

1 > 1 
-Tlx, o ) = r, — | | 



T„,h 


1 

_Lr 


t 





TXx, 0 ) = T, 



T(t) en funcion del Bi 


FlU UK 5.4 Transient tem|>erature distributions for different Biot numbers in a plane wall 
symmetrically cooler! by convection. 

4 4 4 
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Reemplazando L=V/As en la ecuacion de 0 se obtiene 
hAt _ ht _ hL k f _ hLaf 

p V c pc L k pc L 2 k L 2 


hAt 
p V c 


= Bi • Fo 



at 

L 1 


* Fo es el mimero de Fourier. El cual es un tiempo adimensional 
que caracteriza problemas transitorios 
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Bolas de acero al carbon (p=7833kg/m 3 , k=54W/m°C, 

Cp=0.465kJ/kg°C y a=1. 474x1 0 -6 m 2 /s) de 8 mm de diametro se 
recuecen calentandolas primero hasta 900°C en un homo y a 
continuation, dejandolas enfriar con lentitud hasta 100°C en aire 
ambiente a 35°C. Si h=75W/m 2 °C, determine cuanto tardara el proceso 
de recocido. Si se deben recocer 2500 bolas por hora, determine la 
razon total de transferencia de calor de las bolas al aire ambiente. 


Homo 

-,POO n C 

0 

Aifc r 35“C 

IBola dc iffliQ 
f I 

f) o 

|J0 : C 

Q 
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Considere una plancha de 1000 W cuya base esta hecha de aleacion de 
aluminio 2024-T6 (p=2770 Kg/m 3 , Cp=875 J/Kg* °C, a=7.3*10- 5 m 2 /s). La 
placa base tiene un area superficial de 0,03 m 2 y un espesor de 5mm. En un 
principio la placa esta en equilibrio termico con el ambiente a 22°C. Si el 
coeficiente de transferencia de calor por conveccion en la superficie de la 
placa base es de 12 W/m 2 °C y se supone que el 85% del calor generado en 
los alambres de resistencia se transfiere a la placa, determine el tiempo que 
pasara para que la temperatura de esta llegue a 140 °C. <j,Es realista suponer 
que la temperatura de la placa es uniforme en todo momento? 
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• Pared plana T(x,t) 


fi r . e gen _ l dT 


dx- 


+ 


a dt 


Considerar: sin generation, 1 D, prop, fisicas ctes, espesor 
2L, inicialmente a Ti expuesta al ambiente a Too con un h 


d T 1 dT 


dx 2 d. dt 


T(x,0)=T, 




ST(0 ,t) 


dx 


- k — ^ = h[T{L,t ) 


dx 

T—f {x,L,t,k,a,h) 


-T 1 

-*■ ooj 


0, 


awF 


■4 


2L 


Adimensionalizando con X=x/L y 


e=(r-r 00 )/(r,-r 00 ) 
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■ Solucion exacta 

aplicando el metodo de separation de variables con 

Q(X ,Fo)=F (X)-G(Fo) 
imponiendo las condiciones de contorno se obtiene: 

00 2 

e=£c„c“'- F °cos(\„jf) 

n= 1 

4sin \ n 

C " 2\ n +sin(2\„) 

\ n tan \ n = Bi 
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* Solucion aproximada 

si Fo>0.2 la serie se puede aproximar usando solo el primer 
termino (error < 2%) 

0 = C 1 e~ x 1 Fo cos (\X) 

Ci y A1 se obtienen de tablas en funcion del Bi. 

— \ 2 Fo 

* En el centra de la pared 1 0 

* Las temperaturas en cualquier posicion estan relacionadas con la 
del centra asi: 

0 = 0 o cos(A 1 X) 

- La variacion de la temperatura en cualquier punto respecto al 
tiempo es la misma que la del centra de la pared 






Julian Jaramillo 




Escuela Ingenieria Mecanica 
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Solucion aproximada 


para el caso de cilindros y esferas se puede realizar el mismo 
procedimiento. 


* Cilindro: 


®cil 


C 1 e 


Fo 


0 


Jo(\r/r 0 ) , 0 Oci; 
CI -=Jo(\ 1 rlr 0 ) 


C 1 e 


XjFo 


0, cil 


* Esfera: 


^ -X^Fo 


®esf — C x e 


sin (X 1 r/r 0 ) 


\r/r 


0 


9n ^ = 


0 ,esf 


^ —\Fo 

C t e 


Q esf _sin {\r/r 0 ) 


0 


0 ,esf 


\r/r 0 
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Conduccion unidimensional transitorio 
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Solucion aproximada flujo de calor 

Qmax = pVC p (Too- T l) [M] 


Placa plana: 


Q 


placa 


Q 


= 1-0 


sin(A 1 ) 


max , placa 


0 , placa 


\ 


Cilindro: 


Q 


cil 


Q 


= 1-20 




0 , cil 


max, cil 


K 


Esfera: 


Q 


Qesf _ 1 0 sin(\ 1 )-\ 1 cos(\ 1 ) 

- 1 ,esf ,3 


max , esf 


K 
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Ai=Ci 


Fared plana Ciiindro Ester a 


Bi 


A l 

ki 

Aj 


■ 3 j 

O.Ol 

Q .0998 

1 . Q 017 

0.1412 

1.0025 

0.1730 

1.0030 

0.02 

Q . 14 IQ 

1.0033 

0.1995 

1.0050 
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1 .0060 

0.04 
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1.0066 

0.2314 

1.0099 

0.3450 

1 .0120 

0.06 

0.2425 

1.0093 

0.3438 

1.0148 

0.4217 

1.0179 

0.03 

0.2791 

1 . 0 130 

0.3960 

1.0197 

0.4860 

1.0239 

o.i 

0.3111 

1.0161 

0.4417 

1.0246 

0.5423 

1 .0298 

0.2 

0.4328 

1.0311 

0.6170 

1.0483 

0.7593 

1.0592 

0.3 

0.5218 

1.0450 

0.7465 

1.0712 

0.9203 

1 .0380 

0 . 4 - 

0.5932 

1.0580 

0.3516 

1.0931 

1.0528 

1.1164 

0.5 

0.6533 

1.0701 

0.9408 

1.1143 

1.1656 

1.1441 

0.6 

0.7051 

1.0814 

1.0184 

1.1345 

1.2644 

1.1713 

0.7 

0.7506 

1.0918 

1.0873 

1.1539 

1.3525 

1.1978 

0.8 

0.7910 

1.1016 

1.1490 

1.1724 

1.4320 

1 .2236 
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0.3274 

1.1107 

1.2048 

1.1902 

1.5044 

1 .2488 
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1.2558 

1.2071 
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1.0769 
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1.3384 

2.0283 
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3.0 

1.1925 
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1.4191 

2.2889 

1.6227 

4.0 

1.2646 

1.2287 

1.9081 

1.4698 

2.4556 

1.7202 

5.0 

1.3138 

1.2403 

1.9898 

1.5029 

2.5704 

1.7370 

6.0 

1.3496 

1.2479 

2.0490 

1.5253 

2.6537 

1 .3338 

7.0 

1.3766 

1.2532 

2.0937 

1.5411 

2.7165 

1.3673 

a.o 

1.3978 

1.2570 

2.1286 

1.5526 

2.7654 

1 .3920 

9.0 

1.4149 

1.2593 

2.1566 

1.5611 
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1.9106 

10.0 

1.4289 

1.2620 
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1.5677 

2.8363 

1 .9249 

20.0 

1.4961 
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2.2880 

1.5919 

2.9857 

1.9781 

30.0 

1.5202 

1.2717 

2.3261 

1.5973 

3.0372 

1 .9398 

4 - 0.0 

1.5325 

1.2723 

2.3455 

1.5993 

3.0632 

1.9942 

50.0 

1.5400 

1.2727 
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100.0 
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OB 
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1.2732 

2.4048 

1.6021 

3.1416 

2.0000 



Julian Jaramillo 





Escuela Ingenieria Mecanica 


Conduccion unidimensional transitorio 



n 


Jt U> 

0.0 

1.0000 

0-0000 

o . i 
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Funciones de Bessel de cero y 
primer orden 
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• Solucion aproximada 

para Fo>0.2 se desarrollaron graficas que permiten aproximar la 
temperatura y el flujo de calor en placas, cilindros y esferas en 
funcion del Bi y de Fo. 

• Son graficas dificiles de leer para Fo<1 
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Graficas de Heisler 
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To: temperatura en el piano central de la placa=T(x=0,t) 
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= Fo 


T(r, t) — 

Ti-T^ 


"t-t^ 

Jo-T. 


A 

/FIG. 4 


(t -T 
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T -T 

l 1 oo 


\ 

/FIG.3 


To: temperatura en el eje 
del cilindro=T(r=0,t) 
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A 
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To: temperatura en el 
centro de la 

esfera=T(r=0,t) 
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Es un cuerpo idealizado que tiene una sola superficie plana y se 
extiende haci'a el infinito en todas las otras direcciones 

Se usa para indicar que el cambio de temperatura en la region cercana 
a la superficie se debe solo a las condiciones en la superficie plana 

A una cierta profundidad las condiciones en los alrededores no afectan 
la temperatura 

Se considera transferencia de calor unidimensional perpendicular a la 
superficie plana 
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Case (1) 
Tlx, 0) = Ti 
no, t) = t s 



\-*x 


Case (2) 

Tlx, 0) = T t 
-kd17dx\, = 0 = q" o 



h-.t 


Case (3) 

Tlx, 0) = ^ 

-kdT/dx\ x=0 = h[T oo -m, /)] 


T„,h 

n 


I—* 






Figure 5.7 Transient temperature distributions in a semi-infinite solid for three surface 
conditions: constant surface temperature, constant surface heat flux, and surface convection. 
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Aproximacion valida para intervalos de tiempo cortos 
Cambios subitos de las condiciones termicas en la superficie 


8 2 T _ 1 8T 

8x 2 a 8 1 
T(x^co,t) = T i , 


T(0,t) = T s 
T{x,0) = T i 


Haciendo un cambio de variables 


n V4a t 8r\ 2 n dn 

r(n=o)=r s , r(n->oo)=r ( 
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Soluciones 

1 . Temperatur^i superficial conc^cida y (^onstante 


T-T, 


— erfc 


x 


2 Vat 


. q,(t)= 


k(T s -T.) 


tt at 


2. Flujo de calor superficial constante 

Qs 


T{x,t)-T,= 


4a t 


TT 


exp 


x 


4at 


— x erfc 


x 


2 Vat 


3. Conveccion superficial 

T(x,t)-T, ' 

- — erfc 


T -T 

-*-oo -*■ j 


\ 

X 


2 Vat j 



exp 


2 

hx h at 




erfc 


x /Wat 


2 Vat 


\ 




Donde erf es la funcion de error y erfc=1 -erf funcion error complementaria 
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Funcion de error 


w 


erf(w)=—j=j e u du 

vTT o 


Funcion de error complementaria 

erfc{w)= l-erf(w) 


La variable u es una variable ficticia 


X 

erf(x) 

erfc(x) 

X 

erf(x) 

erfc(x) 

0.00 
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0.15 
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0.25 
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1.80 
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0.3286268 

0.6713732 

1.90 

0.9927904 

0.0072096 
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2.00 

0.9953223 

0.0046777 

0.40 

0.4283924 

0.5716076 

2.10 
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0.45 
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0.5245183 

2.20 

0.9981372 

0.0(^8628 

0.50 

0.5204999 

0.4795001 

2.30 

0.9988568 

0.0011432 

0.55 

0.5633234 

0.4366766 

2.40 

0.9993115 

0.0006885 

0.60 

0.6038561 

0.3961439 

2.50 

0.9995930 

0.0004070 

0.65 

0.6420293 

0.3579707 

2.60 

0.9997640 

0.0002360 

0.70 

0.6778012 

0.3221988 

2.70 

0.9998657 

0.0001343 

0.75 

0.7111556 

0.2888444 

2.80 

0.9999250 

0.0000750 

0.80 

0.7421010 

0.2578990 

2.90 

0.9999589 

0.0000411 

0.85 

0.7706681 

0.2293319 

3.00 

0.9999779 

0.0000221 

0.90 

0.7969082 

0.2030918 

3.10 

0.9999884 

0.0000116 

0.95 

0.8208908 

0.1791092 

3.20 

0.9999940 

0.0000060 

1.00 

0.8427008 

0.1572992 

3.30 

0.9999969 

0.0000031 

1.10 

0.8802051 

0.1197949 

3.40 

0.9999985 

0.0000015 

1.20 

0.9103140 

0.0896860 

3.50 

0.9999993 

0.0000007 
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Una tabla gruesa de madera (k=0.17W/m°C y a=1.28 x 
10 -7 m 2 /s) que esta inicialmente a una temperatura uniforme de 
25°C se expone a gases calientes a 550°C durante un periodo 
de 5min. El coeficiente de transferencia de calor entro los gases 
y la tabla es de 35W/m 2 °C. Si la temperatura de ignicion de la 
tabla es de 450°C, determine si se encendera. 
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Conduction en sistemas multidimensionales 
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* Algunas configuraciones pueden ser estudiadas usando los 
procedimientos explicados anteriormente y superponiendolos 
==>solucion producto (multiplicacion de casos ID) 

* <Deben estar rodeados por el mismo fluido a una Too y h constantes 

* Un cilindro corto (2D), por ejemplo, puede ser considerado como una 
interseccion entre un cilindro y una placa plana infinitos que se cruzan 
perpendicularmente 
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* La temperatura en el cilindro corto sera: 


T{r ,x,t) — T 


oo 


T-T 

Z oo 


cil. corto 


T{x,t)-T 

T-T" 


oo 


placa 


T(r,t)-T 

T-T" 


oo 


cil 


* Las soluciones ID se denotan 


QplaJ X >t) = 


T (x,t) — T 


oo 




T-T 


l oo 


placa 


Q J r ’t)= 


T(r,t)-T 


oo 


r,-r 


l oo 


cil 


® semiinf ( ^ 9 ^ ) 


T{x,t)-T 


oo 


T-T 


l oo 


so/, semiinfinito 
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* La transferencia de calor en el cilindro corto sera: 


Q ^ 


Q ^ 

+ 

1 

Q ^ 


1- 

Q ^ 


^ Q max i 

total 2D 

^ Qmax j 

^ Qmax ) 

2 

^ Qmax j 

1 


* Para un cuerpo 3D es: 


Q ' 


\ Q max I 


I ^ \ 


total 2D 


Q 


\ Q max I 


I ^ \ 


Q 


\ Q max I 


1- 


_Q_ 

\ Qmax I 


1 


Q ^ 


1- 

Q ^ 


^ Qmax t 

3 

^ Qmax i 

1 


1- 


_Q_ 

\ Qmax / 
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La temperatura un prisma solido de medidas 2Hx2Lx2W: 


Q prisma [X >y> Z > t) = ® placA* > place { Y , * ) 0 place l Z > f 


<N 



r\... 


0 



h 


H 
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T 


0(r,t) = O cyi (r,t) 
(a) Infinite cylinder 


Lb 


0(r,t)=0 cyl 0 wall (x,t) 
Short Cylinder 



0( r >t) = 0nA e mi-inf(x>t) 

(a) Semi-Infinite cylinder 



Semi-infinite medium 






Quarter-Infinite Medium 


21 


-►T 


e(x,t)=e wall (x,t) 
Infinite Platefor piane Wall) 


Corner Region of a I arge Medium 





0(x.y.t)-0 mM (x.t)O i¥ (y,t) 

Semi-Infinite Plate 
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*V 



Quarter-Infinite Plate 0(x,t )— & S emi-inf( x >^ )@*ull&>0 

Infinite Rectangular bar 



0 w Jx. t )0 wall (y, t (z.t) 

Semi-Infinite Rectangular bar 



O wall (x,t>0 waU ( y.t )O wM ( z.t ) 

Rectangular Parallelepiped 
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Un bloque rectangular de hielo (k=2.22W/m°C y a=0.124xl0 -7 m 2 /s) 
de 5cm de alto y base cuadrada de 4 cm X 4cm inicialmente a -20°C se 
coloca sobre una mesa en una habitacion a 18°C. El coeficiente de 
transferencia de calor sobre las superficies expuestas del bloque de 
hielo es de 12 W/m 2 °C. Si se descarta toda transferencia de calor de la 
base hacia la mesa, determine cuanto tiempo transcurrira antes que el 
hielo se empiece a fundir. <^En donde se formaran las primeras gotas de 
liquido? 
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Considere un bloque cubico cuyos lados tienen 5cm y un bloque 
cilmdrico cuya altura y diametro tambien son de 5cm. Los dos bloques 
se encuentran al principio a 20°C y estan hechos de granito 
(k=2.5W/m°C, alpha=1.15xl0' 6 m2/s). Ahora los dos bloques se 
exponen en un homo a gases calientes a 500°C sobre todas sus 
superficies, con un coeficiente de transferencia de calor de 40W/m2 
°C. Determine la temperatura en el centro de la superficie superior de 
una de las configuraciones al cabo de lOmin. 
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- La disponibilidad de computadores rapidos y software poderosos y 
faciles de usar han tenido un alto impacto en la practica ingenieril 

* Hoy dia los ingenieros tienen que conocer el fenomeno fisico, el 
modelo matematico que lo rige, aplicar correctamente las 
herramientas computacionales y ser capaces de interpretar los 
resultados 

- Es un error usar un codigo sin entender como se esta resolviendo el 
problema ==> el codigo no debe ser una caja negra 

- El know-how adquirido al resolver problemas sencillos es 
indispensable 
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Product 

Specification 



CAD: computer-assisted design 
CAE: computer-assisted engineering 
CSM: computational solid mechanics 
CFD: computational fluid dynamics 
CAA: computational aero-acoustic 
CEM: computational electro-magnetics 
CAM: computer assisted manufacturing 




Efficiency (at cruise) 
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Plate 1.1.4 Impact of CFD on SNECMA fan performance, over a period of 30 
years. From Escuret et al. (1998). 



Engine/airframe integration 

• Simultaneous design 

• Three engine installations 

• Including exhaust effects 
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Chronology and impact 
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Introduction: estructura simulation CFD 
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Real 

World 


Part I 


Levels of Approximation 
and Mathematical 
Model 


Discretization of the 
Mathematical Model 


'« 

■ , it 

Discretization of the 
Space Domain 

Grid Generation _ 

4 

Discretization of the 

Model Equations jl 



Part II 



Part 


Part IV 


Select a discretization 
method and define 
a numerical scheme 


Analyze the 
numerical scheme for 
consistency, stability 
and accuracy 


Solve the resulting 
algebraic system and 
optimize the 
convergence rate 
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* Los metodos numericos convierten las ecuaciones diferenciales en 
ecuaciones algebraicas 


El dominio es dividido en un conjunto de puntos, volumenes, etc. 


* La solucion analitica provee una funcion continua, mientras la solucion 
numerica es discreta 


Plane Wall 
Differential Equation 

d -X + i = o 

dx 2 k 


Valid at every point 
Finite difference equation 


r._,-2r.+r 


m m+ 




/ 



Ax' k 

at discrete point 
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* Los metodos numericos son basicamente de tres tipos: elementos 
finitos, diferencias finitas y volumenes finitos. 

* Elementos finitos generalmente usados en resistencia de materiales y 
diseno de maquinas. 

* En el metodo de las diferencias finitas las derivadas son reemplazadas 
por diferencias. 

* En volumenes finitos se realiza una integracion considerando que el 
valor en el punto es representativo del volumen de control que lo rodea 
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* En diferencias finitas se usa la definition de la derivada 

df(x) = f(x+Ax)-f(x) 
dx Ax 
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